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Abstract — Although  designing  and  tuning  a  proportional- 
integral- derivative  (PID)  controller  appears  to  be  conceptually 
intuitive,  but  it  can  be  hard  in  practice,  if  multiple  (and  often 
conflicting)  objectives  such  as  short  transient  and  high 
stability  are  to  be  achieved.  Traditionally  Ziegler  Nichols  is 
widely  accepted  PID  tuning  method  but  it's  performance  is 
not  accepted  for  systems  where  precise  control  is  required.  To 
overcome  this  problem,  the  online  gain  updating  method 
Augmented  Ziegler-Nichols  PID  (AZNPID)  was  proposed,  with 
the  amelioration  of  Ziegler-Nichols  PID's  (ZNPID's)  tuning 
rule.  This  study  is  further  extension  of  [1]  for  making  the 
scheme  more  generalized.  With  the  help  of  fourth  order 
Runge-Kutta  method,  differential  equations  involved  in  PID 
are  solved  which  significantly  improves  transient  performance 
of  AZNPID  compared  to  ZNPID.  The  proposed  augmented 
ZNPID  (AZNPID)  is  tested  on  various  types  of  linear  processes 
and  shows  improved  performance  over  ZNPID.  The  results  of 
the  proposed  scheme  is  validated  by  simulation  and  also 
verified  experimentally  by  implementing  on  Quanser's  real 
time  servo-based  position  control  system  SRV-02. 

Index  Terms — Discrete  PID  controller,  Ziegler- Nichol's  Auto- 
tuned  PID  controller,  AZNPID 

I.  Introduction 

With  its  three-term  functionality  covering  treatment  to 
both  transient  and  steady-state  responses,  proportional- 
integral-derivative  (PID)  control  offers  the  simplest  and  yet 
most  efficient  solution  to  many  real- world  control  problems. 
Since  the  invention  of  PID  control,  in  1910  and  the  Ziegler- 
Nichols'  (Z-N)  straightforward  tuning  methods  in  1942  [2,  3], 
the  popularity  of  PID  control  has  grown  tremendously.  With 
advances  in  digital  technology,  the  science  of  automatic 
control  now  offers  a  wide  spectrum  of  choices  for  control 
schemes.  However,  an  extensive  survey  on  the  regulatory 
controllers  used  in  industries  reveals  that  97%  of  them  are  of 
PID  structure  [4].  Particularly  at  lowest  level,  as  no  other 
controllers  match  the  simplicity,  clear  functionality, 
applicability,  and  ease  of  use  offered  by  the  PID  controller. 
Its  wide  application  has  stimulated  and  sustained  the 
development  of  various  PID  tuning  techniques,  sophisticated 
software  packages,  and  hardware  modules. 

Though  many  tuning  methods  have  been  proposed  for 
PID  controllers  so  far,  none  of  the  method  could  replace 
Ziegler-Nichols'  (Z-N)  PID  tuning  rule  due  to  its  good  initial 
settings  and  ease  of  use[4].This  rule  performs  satisfactorily 
for  first  order  system,  but  they  fail  to  provide  acceptable 
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performance  for  higher  order  and  non-linear  systems  due  to 
large  overshoots  and  poor  load  regulation  [5]. 

To  overcome  these  drawbacks  many  auto-tuning  methods 
are  proposed  [6-13].  In  [6]  Mixed  Fuzzy-  PID  control  algorithm 
is  applied  for  digital  missile  rudder  servo  system  for  controlling 
two  closed  loops,  position  loop  and  speed  loop.  Paper  [7] 
proposes  a  GPC  (generalized  predictive  control)  implicit  PID 
control  algorithm  applied  to  the  heat  exchanging  station 
control  system,  which  has  features  of  long  time  delay,  time 
varying  and  non-linearity.  Also  in  [8]  PID  is  combined  with 
General  Predictive  Control  (GPC)  such  that  the  computational 
tasks  of  GPC  are  reduced  by  taking  reciprocal  of  number 
instead  of  inverse  of  matrix.  In  [9]  PID  algorithm  is  combined 
with  Model  Predictive  Control  (MPC),  the  restricted  objective 
function  is  altered  to  PID  form.  Through  feed  forward  action 
ameliorative  dynamic  of  control  is  achieved  in  [10]. A 
comparative  study  of  three  PID  algorithms  i.e.  Time 
Constrains,  Observer-based  algorithm  and  Heuristic 
Algorithm  is  done  for  Send  -on-Delta  Sampling  in  [1 1]. Based 
on  the  fractional-order  PID  control  algorithm  and  dynamic 
matrix  control  (DMC)  algorithm,  the  fractional  order  PID 
dynamic  matrix  control  (FOPID-DMC)  algorithm  is  proposed 
in  [12],  for  better  dynamic  performance  and  stability.  In  [13] 
state  of  art  of  PID  control  is  discussed  with  design, 
application,  performance  and  future  scope  of  PID. 

Most  of  the  times  while  working  with  the  real  time  systems, 
dynamics  of  the  system  are  not  completely  known.  In  such 
circumstances  to  deal  with  variable  operating  conditions 
online  gain  modification  scheme  is  proposed  in  [4],  which 
simultaneously  adjusts  proportional,  integral  and  derivative 
gains  depending  on  the  instantaneous  error  (e)  and  change 
of  error  (Ae)  of  the  controlled  variable.  The  study  is  further 
extended  in  [1]  for  getting  better  result's  for  AZNPID 
(Augmented  Ziegler-Nichols  (ZN)  tuned  PID  controllers) 
algorithm  by  application  of  Runge-Kutta  to  solve  differential 
equations  in  PID.  As  mentioned  earlier  this  paper  is  extension 
of  above  study.  In  this  paper  modified  AZNPID  algorithm  is 
applied  to  second,  third  and  fourth  order  processes  along 
with  its  hardware  verification.  Stability  aspects  of  processes 
are  also  discussed  with  the  help  of  Nyquist  stability  criteria. 

The  rest  of  the  paper  is  organized  as:  Section  II  gives  an 
overview  of  PID  and  Auto-tuned  PID  controller.  Section  III, 
illustrates  the  experimental  results  for  different  systems.  In 
Section  IV,  results  for  application  of  AZNPID  on  QUANSER's 
SRV-02  plant  are  discussed. 


16 


^cAC 


Full  Paper 


ACEEE  Int.  J.  on  Control  System  and  Instrumentation,  Vol.  4,  No.  1,  Feb  2013 


II.  Designing  of  Auto-Tuned  PID  Controller 

The  Conventional  PID  is  shown  in  Fig.  1. For  hardware 
implementation,  discrete  PID  controller  is  considered. 

A.  Discrete  PID  Controller 

Discrete  form  of  a  conventional  ZNPID  can  be  described 
as, 

u(k)  =  KP  [  e(k)  +  -  e(j)  +  -  [e(k)  -  e(k  - 1)  ]]         (i) 

ti  ts 

Where,  ts  is  sampling  period,  n  is  index  and  j  is  referring  to 
time  constant.  KP  is  proportional  gain,  ti  is  integral  time 

constant  or  reset  time,  td  is  derivative  is  derivative  time 
constant  and  e  is  the  error  between  the  reference  and  process 
output  y  for  instance  k.  Discrete  time  integral  gain  J£i  and 
derivative  gain  Kd  obtained  as, 

Ka  =  KP  —  and  Kd  =  Kp  — 

ti  ts 

The  values  of  Kp ,  KandKd  are  decided  using  any 
standard  rule.  For  this  study,  we  have  used  Ziegler-Nichol's 
Tuning  rule  [14].  The  values  of  critical  gain  Kcr  and 
corresponding  time  period  pcr  are  determined  experimentally 
or  can  be  calculated  using  bode  plot.  Form  those 
values  KP  ,ti  and  td  will  be  determined  as  given  in  Table  I. 

Table  I.  Ziegler-Nichol  Tuning  Rule 


Gain  Coefficient 

PID 

Proportional  Gain  (  KP  ) 

0.6  Kcr 

Integral  time  ( ti ) 

0.5  Per 

Derivative  time  ( td  ) 

0.125  Pcr 
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Figure  1.  Block  Diagram  of  PID 
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Figure  2.  Block  Diagram  of  AZNPID 


B.  Augmented  Ziegler-Nichols  PID  Controller 

Fig.  2  describes  the  basic  structure  of  AZNPID.  We  get 
initial  values  of  KP ,  K anc*  Kd  from  Table  I.  The  e(k)  and 
Ae(k)  is  calculated  as  [4] 

e(k)  =  ref-y(k) 

e(k) 


eN(k)  = 


-) 


(2) 
(3) 


abs(refn) 
AeN(k)=eN(k)-eN(k-l) 
P(k)  =  eN(k)*AeN(k) 
Values  of  KP  ,  Ki  anc*  Kd  are  modified  depending  on  value 
of  f3(k)  after  every  sample  instant  as, 

Kpm=Kp(l  + *,.!£(*)  I) 
K1-=K(l  +  *u^(ik)l) 

K,m  =  m\+k*\p(k)\) 

k 


(4) 
(5) 


(6) 
(7) 

(8) 


u-o  (k)  =KPm  (k)e(k)+Kim  (k)£e(i) 

(i=0) 

+  Kdm(k)Ae(k) 


(9) 


where, 


Ae(k)  =  e(k)-e(k-l) 

e(j)  =  e(k)-e(k-l) 

In  (6)  to  (9)  kH  ,  ki2 '  kn  are  three  positive  constants  which 
gives  variation  in  process  value  to  achieve  desired  output. 
For  this  study,  performance  indices  are  chosen  as  kn  =1, 
ki2  =1»  ku  =12.Equation  (2)  to  (5)  indicates,  when  process 

output    is  less  than  set  point  (refn),  value  of  KPmis  is 

greater  than  KP  so  that  the  process  achieves  set-point 

quickly,  Kd  m  is  greater  than  Kd  so  as  to  avoid  oscillations  in 

process  output.  At  the  same  time  j^i m  is  less  man  K t0  cor" 
rect  offset  errors,  which  plays  a  minor  role  when  the  process 
is  reaching  to  set  point.  As  the  process  reaches  to  set-point 

values   ofKPm     Km'Kdmare   approximately   equal 

to  KP  ,  K  » Kd  -As  the  process  output  increases  beyond  set 

point,  the  values  of,  KP  m  ,  K  m  » Kd m  increases,  this  eventu- 
ally increases  control  u  ( Uauto  >  u ).  This  brings  back  pro- 
cess to  desired  set-point. 

A.  Steps  for  implementation  of  AZNPID  algorithm: 

Stepl:  Obtain  gain  margin  and  phase  margin  of  system 
either  using  bode  plot  method  or  Ziegler  Nichol's  -Continuous 
Oscillation  Method. 
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Step2:  In  ZN  Continuous  Oscillation  Method,  critical  gain  of 
system  is  Kcr  and  its  corresponding  period  of  oscillation  is 

Pcr 

Step3 :  Applying  Zeigler-Nichols  formula,  the  values  of 

KPm  Kim'Kdm  is  obtained  from  table  I 
Step4:      Choosing      the      appropriate      values      of 
constants  kn  ,  k12  and  AZNPID  is  applied  to  the  system  using 
(6)  to  (9).  K13 

III.  Experimental  Results 

The  robustness  and  stability  of  AZNPID  is  tested  for  second 
order,  third  order  and  fourth  order  systems. 

A.  Second  Order  System 

The  transfer  function  for  second  order  system,  with  dead  - 
time  is  given  as 


-0.2s 


Gp  = 


/    /1 ^  (10) 

(s(l  +  s)) 

Referring  Fig.  3,  one  can  observe  that,  after  application  of 
AZNPID  there  is  considerable  reduction  in  peak  overshoot 
of  the  second  order  system.  For  testing  purpose,  the  noise 
with  0.1%  magnitude  and  sinusoidal  disturbance  of  amplitude 
1  with  frequency  1Hz  is  applied  to  system  and  response  is 
observed  in  Fig.  4.  Stability  of  the  system  is  analyzed  by 
using  Nyquist  plot,  for  the  above  system  is  shown  in  Fig.  5. 

B.  Third  Order  System 

Consider  a  third  order  non-minimum  process, 


GP  = 


(1-O.ls) 


(11) 


a+sr 

From  Fig.  6,  one  can  observe  that  there  is  considerable 
improvement  in  the  transient  performance  of  the  system  with 
AZNPID  compared  to  ZNPID.  Fig.  7  shows  effect  of  noise 
and  external  disturbances  for  third  order  system.  Fig.  8  shows 
stability  analysis  of  the  system  using  Nyquist  plot.  Fig.  12 
shows  effect  of  variation  of  performance  indices  for  the 
system  above. 

C.  Fourth  Order  System 

Consider  the  fourth  order  linear  system  as 

(7s +  1) 


G  P  = 


(s4+10s3+35s2+50s  +  24) 


(12) 


For  fourth  order  system,  from  Fig.  9,  one  can  observe 
that,  as  order  of  system  increases,  oscillations  of  ZNPID  also 
increased  but  there  is  significant  effect  of  modified  derivative 

gain  Kd m  m  AZNPID.  Due  to  this  there  is  also  a  considerable 
change  in  settling  time  and  peak  overshoot.  Fig.  10  shows 
effect  of  noise  and  external  disturbances  for  fourth  order 
system.  Fig.  1 1  shows  stability  analysis  of  the  system  using 
Nyquist  plot. 
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Figure  3.  Comparison  of  responses  of  PID  and  AZNPID  for  Second 
order  system 


15  I        " " ■ 1 1 


'• ! 


ZNPKJ 

AZNPID 


10  15  20 


:r- 


30 


Figure  4.  Comparison  of  responses  of  PID  and  AZNPID  in  presence 
of  a)  noise  and  b)  disturbance  for  Second  order  system 

Nyquifi  Diagram 


Figure  5.  Nyquist  Diagram  for  Second  Order  System 

IV.  Case  Study 

To  demonstrate  the  performance  of  the  proposed 
algorithm,  Quanser 's  SRV-02  plant  is  used  for  experimentation. 
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Figure  6.  Comparison  of  responses  of  PID  and  AZNPID  for  Third 
order  system 
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Figure  7.  Comparison  of  responses  of  PID  and  AZNPID  in  presence 
of  a)  noise  and  b)disturbance  for  Third  order  system 


Figure  8.  Nyquist  Diagram  for  Third  Order  System 

The  digital  servo  rig  has  two  parts,  namely  hardware  and 
software.  The  hardware  units  are  mechanical  units  and  digital 
unit.  Mechanical  unit  consists  of  power  amplifier,  dc  servo 
motor  and  tacho-generator,  absolute  and  incremental  digital 
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Figure  9.  Comparison  of  responses  of  PID  and  AZNPID  for  Fourth 
order  system 
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Figure  10.  Comparison  of  responses  of  PID  and  AZNPID  in 
presence  of  a)  noise  and  b)  disturbance  for  Fourth  order  system 
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Figure  11.  Nyquist  Diagram  for  Fourth  Order  System 

encoders,  input  and  output  potentiometers  [15,  16].  The 
digital  unit  carries  ADC  and  DAC  for  signal  conversion, 
switching,  multiplexing,  encoder  circuits,  and  PWM  motor 
drive.  With  the  help  of  SIMULINK  and  Embedded  MATLAB 
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Figure  12.  Comparison  of  responses  of  AZNPID  for  different 
values  of  performance  indices  for  Third  Order  system 

function,  under  MATLAB  7.3  (2009b)  control  algorithm  was 
implemented. 

A.  Mathematical  Model 

In  the  Quanser's  servo  plant  module,  servo  motor  is 
connected  to  optical  encoder  through  gearbox.  Gear  ratio 
can  be  changed  according  to  requirement  of  user  [15, 
16]  .Transfer  function  of  servo  plant  can  be  written  as, 


0(s) 


K 


(14) 


V„,(S)  S(ST+1) 

Where,  0(s)  is  Laplace  transform  of  position  of  the  load 

shaft,  Vm(s)  is  Laplace  transform  of  input  voltage,  K  is  the 

steady-state  gain,  %  is  the  time  constant,  s  is  the  Laplace 
operator. 

State  Space  representation  of  Servo  Plant  was  obtained 
with  position  and  velocity  as  states  as, 

"0 
1 


0(s) 


0in(s) 


1 


6{s) 

COm(s) 


+ 


0 

Ky 


Vm(s) 


(15) 


Figure  13.  Experimental  Setup  for  SRV-02 
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Where,  &)m(s)  is  Laplace  transform  of  speed  of  load  shaft. 

By  substituting  specifications  of  the  servo  plant  [15,  16] 
mathematical  model  was  obtained  as, 

210 


GP  = 


s(s  +  10) 


B.  Results 

For  experimentation,  reference  square  wave  was  taken 
from  SIMULINK.  The  experimental  setup  is  shown  in  Fig.  13. 
Here,  the  proposed  controller  was  implemented  using 
SIMULINK  programming.  With  the  help  of  Real  Time 
Workshop  (RTW)  and  Quarc4  Windows  Target  (RTWT) 
environment  the  performance  of  the  proposed  auto-tuning 
scheme  has  been  implemented  for  achieving  the  desired 
position  (given  by  the  reference  signal)  by  the  dc  servo  motor. 

RTWT  communicates  with  the  control  program,  and 
interfaces  the  mechanical  unit  through  the  digital  board.  RTW 
generates  C  code  using  Microsoft  C+  +  Professional  from 
the  SIMULNK  block  diagram  and  acts  as  the  intermediary  for 
two  way  data  flow  from  the  physical  servo  system  to  and  from 
the  SIMULNK  model.  The  tuning  of  the  PID  controller  is  done 
according  to  ZN  relation  (proportional,  integral,  and  derivative 
parameters  are  given  by  the  manufacturer  feedback). 

Fig.  14,  Fig.  15  and  Fig.  16,  show  the  responses  and 
variation  of  control  signals  for  ZNPID  and  AZNPID, 
respectively  for  both  constant  set-point  and  variable  set- 
point  (square  wave). 

One  can  observe  that  the  proposed  scheme  supposed  to 
provide  significantly  improved  performances  for  those 
systems  which  produce  large  overshoots  under  ZNPID. 
However,  the  real  servo  system  does  not  exhibit  such 
behavior;  as  a  result,  we  find  marginal  improvement  in  the 
performance  of  AZNPID  over  ZNPID  as  shown  in  Fig.  14, 
Fig.  15  and  Fig.  16. 

Fig.  17  represents  output  of  AZNPID  in  presence  of  noise 
and  disturbance.  In  sub-plot  2, one  can  observe  chattering 
state  due  to  noise.  In  sub-plot  3, sine  disturbances  are  applied 
to  AZNPID.  In  general,  PID  doesn't  withstand  in  presence  of 
sinusoidal  disturbances,  although  the  plot  is  showing  same 
property  but  the  amplitude  of  diversion  from  reference  is 
very  less. 
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Figure  14.    Responses  of  PID  and  AZNPID  for  constant  reference 
sources  a)  Process  Output  b)  Control  Input 
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Table  II.  Comparison  Of  Responses  For  PID  And  AZNPID   In  Matlab 


Systim 

PBrfonnaiE? 
Psr-sia^tsrs 

ZNPID 

AZNPID 

fl-^ 

%  overshoot 

50. 8 

34.44 

Ris^Tim^EK) 

0.25 

0.472 

^"stl  +  s) 

Sa-tilinE  Time(5ec) 

6.5 

5.2 

%  ov=r shoot 

3,6.5 

24.9 

Ris  =  Tim2-{s5cJ) 

1  416 

1221 

S  a-ttLiiiE  Tim*£s  ac) 

13.52 

7.053 

G_P 

+  3SffT+5fe+24) 

%  overshoot 

56.2 

1109 

RisaTima^se*) 

0.327 

0.346 

S  3-ttliiiE  Tima{s  k) 

15 

6.E49 
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Figure  15.  Responses  of  PID  and  AZNPID  for  square  wave  as 
reference  source:  nrocess  outout 
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Figure  16.  Responses  of  PID  and  AZNPID  for  square  wave  as 
reference  source:  control  input 

Conclusions 

In  this  paper,  a  simple  model-independent,  discrete  auto- 
tuned  scheme  of  PID  controller  is  presented.  Potency  of  the 
algorithm  is  verified  by  applying  on  second  order  process 
with  dead  time,  third  order  non-minimum  process  and  fourth 
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Figure  17.  Comparison  responses  of  AZNPID  a)  response  of 
AZNPID  b)  in  case  of  noise  and  c)  in  case  of  sine  disturbance 

order  linear  process.  Using  proposed  AZNPID  algorithm,  it 
is  observed  that,  the  transient  performance  of  the  system 
improved  significantly,  even  if,  there  is  change  in  the  order  of 
the  system.  The  scheme  adjusts  gain  of  P,  I  and  D  after  every 
sample  instant  based  on  gain  updating  factor  beta  (P). 
Although  variation  of  beta  is  non-linear,  original  structure  of 
PID  remains  intact.  Robustness  of  AZNPID  and  ZNPID  was 
tested  by  adding  noise  and  external  disturbances.  Though 
AZNPID  is  able  to  improve  the  transient  performance,  steady 
state  performance  remains  unaffected.  Thus,  AZNPID 
controller  can  be  used  where  the  system  has  poor  transient 
response  (non-minimum  phase  system,  higher  order  systems). 
To  test  the  efficacy  of  hardware  implementation,  the 
algorithm  is  applied  to  Quanser's  DC-  servomotor  position 
control  system  (SEV-02).  AZNPID  intensifies  the  results  in 
different  circumstances.  Moreover,  the  scheme  can  be  apt 
for  embedding  on  dedicated  hardware. 
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